Manifestation of quantum chaos on scattering techniques: application to low-energy 
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Intensities of LEED and PED are analyzed from a statistical point of view. The probability 
distribution is compared with a Porter-Thomas law, characteristic of a chaotic quantum system. 
The agreement obtained is understood in terms of analogies between simple models and Berry's 
conjecture for a typical wavefunction of a chaotic system. The consequences of this behaviour 
on surface structural analysis are qualitatively discussed by looking at the behaviour of standard 
correlation factors. 
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There is a continuous interest in understanding the re- 
lationship between chaos and quantum mechanics. Long 
time ago, Wigner investigated the influence of chaos on 
quantum, mechanical scattering experiments in nuclear 
systems!!! Thereafter, much work has concentrated on 
the analysis of energy levels of bound states inside closed 
systems (like various types of billiard geometries) . While 
these studies offer obvious advantages, a great deal of 
information is lost by neglecting the examination of the 
wavefunctions. In fact, a good understanding of wave- 
functions is crucial to explain open systems, like the stan- 
dard probe-target-detector setup used in most scattering 
experiments. Therefore, it is quite perplexing to find so 
few examples in the literature related to quantum chaos 
manifested in experiments where wavefunctions are ana- 
lyzed, which should, be emphasized, only correspond to 
closed geometries. tffi In this work, we show that reflected 
intensities of surface scattering experiments, which are 
directly related to the modulus squared of the wavefunc- 
tion, are consistent with quantum chaos. Therefore, we 
are proposing a new class of simple experimental systems 
where quantum chaos is manifested in the properties of 
wavefunctions. 

One reason to expect quantum chaotic behaviour in a 
scattering experiment comes from the existence of clas- 
sical chacp when three or more scattering potentials are 
involved.D Actually, Mucciolo et alJj have recently shown 
that the high energy region of the calculated band struc- 
ture of Si and Al x Gai_ x As is complex enough to obey the 
statistical distribution gf levels corresponding to random 
matrix theory (RMT).Q Based on this statistical analy- 
sis, the authors claim that these systems exhibit quan- 
tum chaos. Although just a theoretical prediction, this is 
a remarkable result because no disorder or incommensu- 
rate geometries are involved, and the physical reason for 
chaos should be found elsewhere (e.g., the intrinsic mul- 
tiple scattering (MS) originating Bloch states and giving 
rise to the crystal band structure). 

In this letter, the manifestation of chaos on standard 
surface scattering techniques like Low Energy Electron 
Diffraction (LEED) or Photoelectron Diffraction (PED) 
is investigated. These techniques are dominated in most 
experimental systems by MS, yielding a clear similitude 
to the band structure problem mentioned above. In sur- 
face structural work, it is a common belief that MS intro- 
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duces a richier but more difficult analysis. This general 
statement is analyzed here from the point of view of the 
chaotic component of the LEED and PED experiments, 
while we draw new attention to a class of quantum syi 
terns obeying a Porter-Thomas probability distribution 
Ultimately, our aim is the understanding of the rela- 
tionship between complex scattering phenomena and the 
emergence of quantum chaos. 

To characterize chaotic wavefunctions Porter and 
ThomasQ advanced the hypothesis that wavefunctions of 
a chaotic system should display a \l statistical probabil- 
ity distribution. Subsequently, this hypothesis has been, 
rigorously justified using the supersymmetry formalism,!!! 
and has been used as a convenient definition of quantum 
chaos, that at-.least can be thought as a necessary con- 
dition. DysonO demonstrated that within the RMT only 
three universal classes can exist depending on whether 
the hamiltonian is constructed with real numbers, com- 
plex numbers or quaternions, corresponding respectively 
to v = 1,2, and 4 degrees of freedom. Since scattering 
wavefunctions are complex numbers, the statistics corre- 
sponding to the universality class v = 2 is expected. 

An interesting theoretical result onJJae wavefunction 
of a chaotic system is due to Berry.E3 Analysing the 
semiclassical mechanics of regular and irregular motion, 
he realized that a typical chaotic wavefunction should 
be a linear combination of plane waves with random k- 
orientations, (at a fixed constant energy), and random 
complex coefficients: 
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Heller et al.0 have further investigated the properties 
of these chaotic wavefunctions, finding that in 2D they 
present characteristic scars. Berry's chaotic wavefunc- 
tion can be interpreted as the result of multiple random 
reflections of a plane wave. It can also be thought as 
the superposition of plane waves originating at N points 
propagating with the same energy in random orienta- 
tions, and mixed with appropriated coefficients. Guided 
by these images, we try to find a physical system where 
a similar wavefunction can be realized. 

First of all, we consider a PED experiment where an 
electron inside an atomic core is excited by an incident X- 
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ray photon, to be subsequently diffracted by a cluster of 
n atoms surrounding the original source. This is a com- 
plicated scenario, but making some approximations that 
have been tested in the literature, a simple-expression for 
the wavefield at a distance R (far-field) is:liil 
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where the complex coefficients, c? u , include appropriated 
scattering factors and the expansion can be extended to 
any desired order of scattering. It is important to realize 
the similarity of this MS series with Berry's one: taking 
away the prefactor and the source wave, and given a fixed 
direction in real space determined by k, it is written as 
a sum of N plane waves e lk ' Tu ' with complex coefficients, 
where the many r u may result oriented in uncorrelated 
directions if enough scattering is allowed. The following 
question arises in this context: how many plane waves 
are necessary to allow Eq. (1) to follow a Porter-Thomas 
statistical distribution? As an example, if k = 2-7T, we 
find that just N = 10 are enough to find a reasonable 
agreement (e.g., see Fig. 1); the result for as few as 
N = 2 is given in the same figure for comparison. 

Secondly, we notice that an expression that is formally 
similar to Eq. (2) can be written for the Diffuse LEED 
(DLEED) wavefield:E3 
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where F a represent generalized scattering factors includ- 
ing MS contributions. Standard LEED I(V) curves can 
be described in the same way, for a given energy, just 
keeping in mind that if the system exhibits perfect pe- 
riodicity in the parallel direction, only a discrete set of 
points given by Bragg conditions are available. 

Before trying to analyze real experiments, a set of con- 
trolled theoretical simulations of relevant systems is con- 
sidered. We investigate the behaviour of the single scat- 
tering term in Eq. (2) performing the summation over 
a set of 500 atoms randomly distributed around the ori- 
gin between r a — 10 a.u. and r a = 150 a.u. (atomic 
units will be used throughout the paper, expressing dis- 
tances in Bohrs and energies in Hartrees). The central 
region (magnified ten times) of a typical | ipk | 2 {k = 6 
a.u.) measured on a sphere at an asymptotic distance is 
shown in Fig. (2). A typical worm-like image is obtained 
when the pattern saturates at high energies or high r a 
distances. The corresponding probability distribution of 
intensities (normalized to the average) is seen to follow a 
X2 distribution rather well. 

Using the same model, we explore the probability dis- 
tribution of intensities produced by a small quasi-regular 
cluster of atoms. A cube of 2 x 2 x 2 Ni atoms centered 
around the origin at a distance of 5 a.u is chosen. Phase 
shifts up to l,nax = 7 are used to compute the scattering 



factors. These phase shifts have been obtained from a 
muffin-tin model, and can be used to represent the scat- 
tering of electrons by atoms in the real structure. In or- 
der to simulate small geometrical irregularities caused by 
relaxations, reconstructions or simply the effect of tem- 
perature, the atoms are randomly displaced from their 
perfect positions in the cube with values uniformly dis- 
tributed between and \/3d. Fig. 1 shows the result of 
such a simulation for d — 1 a.u.. The agreement with the 
Porter-Thomas law is really good. In the same figure, the 
effect of a smaller arbitrary displacement (d = 0.2 a.u.) 
is also shown, together with a similar calculation for a 
double-scattering term in a cubic array of 3 x 3 x 3 Ni 
atoms with comparable results. The analysis of Eq. (1) 
proves that fluctuations are responsible for the appear- 
ance of the ideal Porter-Thomas distribution. Therefore, 
the deviations from the x| distribution observed for the 
scattering series should be explained by studying their 
fluctuations, which is beyond the scope of this work. _ 

We compute full dynamical diffuse Leed intensities!!^ 
for a realistic adsorption geometry on the system 
O/Ni(100) (oxygen is placed on the hollow site at 1.51 
a.u. from the surface). All parameters are takep-,from 
a detailed structural analysis of the same system.O It is 
worthwhile to notice two points: i) All the Ni atoms are 
placed at ideal bulk-like positions. Therefore, there is 
not geometrical disorder in the problem, the main source 
of complexity being MS by the atoms in the ordered lat- 
tice; and, ii) All the calculations are performed at T = 
K, although attenuation effects are taken into account in 
this formalism via an imaginary part (Vbj = 0.15 a.u.) 
added to the energy. Fig. 3 shows the probability dis- 
tribution for intensities calculated theoretically at three 
different energies (12, 14 and 16 a.u.), yielding a similar 
agreement to the Porter- Thomas probability distribution 
as the previous PED example. The same results are ex- 
pected from the analysis of experimental intensities. As 
an example, Fig. 3 includes a single energy (11.1 a.u.) 
extracted from the jexperimental database measured by 
the Erlangen group £j 

The same Porter- Thomas probability distribution 
should also appear when conventional LEED I(V) curves 
are analyzed, because our arguments above are valid 
for any energy. We have simulated theoretically the 
LEED I(V) curvesH for three materials with very dif- 
ferent structures: Cu(100), W(100), and Si(lll). An 
arbitrary non- normal incidence angle (9 = 20°, cj> = 30°) 
breaking the symmetry is chosen. This yields the max- 
imum number of inequivalent beams increasing the sta- 
tistical confidence of the analysis. The first 9 emergent 
beams are used for Cu and W while the first 13 were 
chosen for Si. The energy ranged from 50 to 450 eV for 
both metals and from 30 to 300 eV for the semiconductor, 
yielding a database of about 100 a.u., which is enough for 
our purposes and easily accessible to experiment. The 
imaginary part of the energy is fixed to a constant value 
of 0.15 a.u., and T = K is used again. Finally, we ana- 
lyze the experimental database for c(8 x 2) GaAs(100),E^I 



2 



formed by 13 different beams measured at normal inci- 
dence. Our results are shown in Fig. 4, displaying an 
agreement with the Porter-Thomas probability distribu- 
tion similar to the other examples. 

Guided by these results, we predict the existence of a 
region (II) in parameter space, P, where small changes 
in p (each component defining a relevant parameter for 
the structure), result in rapid changes of the wavefunc- 
tion. On intuitive grounds, it can be assumed that these 
changes must separate exponentially. This region is in- 
termediate between two different ones: a perturbative 
region (I), as for sufficiently small changes in p we expect 
perturbation theory to give a reasonable answer,E3 and 
a random region (III) where wavefunctions for different 
structures are uncorrelated. The existence of these three 
regions is checked by analyzing two R-factors commonly 
used in surface,structure analysis: (i) A root mean square 
displacement!!^ adapted for DLEED (R 2 ) and, (ii) the 
Pendry R-factorE3 (Rp) often used with standard LEED. 

We apply R 2 to compare theoretically calculated 
DLEED intensities for O/Ni(100) as a function of the 
adsorption height h. Arbitrarily, h = 1.51 a.u. is chosen 
as the reference. This is a clean and controlled theoret- 
ical experiment where only the position of one atom is 
changed, but to stay closer to reality, we also consider 
the behaviour of the Rp in a recent structural analysis of 
c(2 x 2)Si/Cu(110) I(V) curves, where, the relaxation of 
the whole surface layer is considered. EJ Both cases show 
the existence of the three regions mentioned above. Re- 
gion I is obviously well characterized by the existence of 
a minimum that imposes a quadratic dependence. Re- 
gion III is also easily identified by the saturation of the 
R- factor: for Rp this happens by construction around 
1 (the maximum value for Rp can be twice this value, 
but saturation starts at values greater than 0.6). For R 2 
we observe that saturation occurs around the value ob- 
tained by comparing two sets of N random intensities, 
so the R-factor is normalized to this value. Region II 
may be characterized by plotting the ln(i?) versus a rel- 
evant component of p, and identifying the interval where 
it behaves like a straight line. When we use these ideas 
to analyze the data, we find that the quadratic region 
(I) extends approximately 0.04 A for the DLEED case 
(i? 2 < 0.2) and 0.05 A for the conventional LEED analy- 
sis (Rp < 0.25). The exponential region (II) extends also 
0.04 A for DLEED (R 2 < 0.5), while it goes to 0.09 A for 
the LEED experiment (Rp < 0.6). Finally, an uncorre- 
lated region extends beyond these intervals, provided we 
do not approach a multiple coincidence minima. We no- 
tice that a perturbative technique where the perturbation 
in the potential is proportional to the atomic displace- 
ments (the so-called Tensor LEED first approximation!!!!) 
is known to break down beyond w 0.1 A. This is close to 
regions I and II, considered respectively a truly pertur- 
bative region (quadratic) and the onset of the breakdown 
for the perturbative approach (exponential). These find- 
ings should bring more rigour to the standard R-factor 
analysis because our analysis allow to identify regions II 



and III, where correlations must be taken as spurious, 
and offer a new theoretical explanation for the definitive 
failure of simple scattering methods in LEED. 

We have analyzed typical wavefunctions for LEED and 
PED experiments in the light of Berry's proposal for 
a generic chaotic wavefunction. Our statistical analy- 
sis shows that scattering wavefunctions computed from 
several models (including perfectly ordered structures) 
follow the Porter-Thomas xi distribution. This property 
is also obtained analyzing experimental data for LEED 
and DLEED. The physical origin of this behaviour is the 
complexity of the scattering. Attenuation effects taken 
into account via a complex optical potential fitted to ex- 
periments, and defects (relaxations or reconstructions) 
do not change this conclusion. Finally, analysing the be- 
haviour of two different R-factors, we have argued the 
existence of three distinct regions, showing the rationale 
behind widely used rules about which R-factors are ac- 
ceptable in standard structural work and which are not. 

We are grateful to Prof. K. Heinz for making available 
to us his experimental DLEED data on O/Ni(100). This 
work has been supported by the CICYT under contracts 
num. PB94-053 and MAT94-0058. 
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FIG. 1. Probability distribution for simplified models: (a) 
Thick solid line: xl l aw i (b) Thin solid line: Berry's wavefunc- 
tion with N — 2; (c) dashed-dotted line: Berry's wavefunction 
with N = 10; (d) dotted line: simple PED model for eight 
Ni atoms on a cube, simple scattering (E — 1.8 a.u., d = 1.0 
a.u.); (e) dashed line: same as (d) with d = 0.2 a.u.; (f) 
dashed two-dotted line: simple PED model for 27 Ni atoms, 
double scattering (E = 1.8 a.u., d = 1.0 a.u.). 




FIG. 2. Small portion of a wavefunction (modulus 
squared) obtained using Eq. (2) for a set of 500 atoms at 
random positions. 
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FIG. 3. Probability distribution for MS intensities in a 
DLEED model. Different energies are shown, (a) Thick solid 
line: xl law; (b) dotted: E — 12 a.u.; (c) dashed: E — 14 
a.u.; (d) dashed dotted: E — 16 a.u.; (e) dashed two-dotted: 
E — 11.1 a.u. (experimental). 
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Figure 4 
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FIG. 4. Probability distribution for MS intensities in a 
LEED model. Different materials are shown, (a) Thick solid 
line: Porter- Thomas law; (b) dotted: Cu(100); (c) dashed: 
W(100); (d) dashed dotted: Si(lll); (e) dashed two-dotted: 
c(8 x 2)GaAs(100) (experimental). 
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